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Abstract 

We investigate dominant late-time tail behaviors of massive scalar fields in 
nearly extreme Reissner-Nordstrom background. It is shown that the oscilla- 
tory tail of the scalar fields has the decay rate of t -5 / 6 at asymptotically late 
times. The physical mechanism by which the asymptotic i~ 5 / 6 tail is gener- 
ated and the relation between the field mass and the time scale when the tail 
begins to dominate, are discussed in terms of resonance backscattering due to 

spacetime curvature. 
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I. INTRODUCTION 



Various interactions of black holes with scalar fields have been extensively studied for 
a long time. Though previous works have been mainly concerned with the evolution of a 
massless scalar field, the analysis of massive ones will be also physically important. For 
example, in higher- dimensional theories, the Fourier modes of a massless scalar field behave 
like massive fields known as Kaluza-Klein modes, and the recent development of the Kaluza- 
Klein idea (e.g., the Randall-Sundrum model [|I[] in string theory) strongly motivates us to 
understand the evolutional features due to the field mass in detail. 

Massive scalar fields in black hole spacetimes can cause interesting phenomena which are 
qualitatively different from massless case. The remarkable example is the vacuum polariza- 
tion ((f) 2 ) of a quantum massive scalar field <p, which is in thermal equilibrium with a nearly 
extreme Reissner-Nordstrom black hole || . Because the black-hole temperature is very low, 
the mass-induced excitation of (0 2 ) results in the resonant amplification at mM ~ 0(1) for 
the field mass m and the black-hole mass M. 

Such a resonance behavior due to the mass of a field interacting with a black hole may 
appear in various processes as a basic feature of black-hole geometry. For a step to support 
this conjecture, in this paper, we turn our attention to a problem on time evolution of 
classical fields. 

The evolution of a massive scalar field in Schwarzschild background was first analyzed 
by Starobinskii and Novikov, using the complex plane approach fU, and they found that 
because of the mass term, there are poles in the complex plane closer to the real axis than 
in the massless case, which lead to inverse power law behavior with smaller indices than 
massless case. Recently, it was pointed out that the late-time tails of massive scalar fields 
in Reissner-Nordstrom spacetime are quite different from massless fields in the existence of 
the intermediate late time tails If the field mass m is small, namely mM -C 1, the 
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oscillatory inverse power-law behavior 

ip ~ sin(mt), (1) 

dominates as the intermediate late-time tails. Note that massless fields decay more rapidly 
without any oscillation, as was first studied by Price ]E[]. The analytical approximation ([!]) 
was derived from the flat-spacetime approximation, where the effects of spacetime curvature 
are neglected. Though the behavior ([!]) was numerically verified at intermediate late times, 
Hod and Piran [||] also suggested that another wave pattern dominates at very late times, 
namely the intermediate tails (0) are not a final asymptotic behaviors, and mentioned that 
"fields decay at late times slower than any power law" . Though similar numerical results in 
the Schwarzschild case were reported by Burko ||, the evolution of massive scalar fields at 
asymptotic late times has been claimed to be inverse power-law behavior. 

Our purpose of this paper is to clarify what kind of mass-induced behaviors dominates 
in the asymptotic late-time tails as a result of interaction of massive scalar fields with a 
black hole. In Sec. [TT] we introduce the black-hole Green's function using the spectral 
decomposition method |7| . (Another new method to treat late-time tails which is called "late 
time expansions" was also proposed recently ||.) In Sec. |TTT| we consider a nearly extreme 
Reissner-Nordstrom background, motivated by the fact that the resonance phenomena in the 
vacuum polarization of quantum scalar fields become clear. Then, based on the procedure 
of asymptotic matching, we construct approximate solutions in the nearly extreme limit. 
In Sec. [IV] we study the intermediate tail which appears in the case of small mass field, 
confirming our result in comparison with (|l|). In Sec. [V| we derive the asymptotic late-time 
tail which is the main result in this paper. The final section is devoted to summary and 
discussion. 



II. GREEN'S FUNCTION ANALYSIS 
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A. Massive scalar fields in Reissner-Nordstrom geometry 

We consider time evolution of a massive scalar field in Reissner-Nordstrom background 
with mass M and charge Q. The metric is 



ds = ~ [ 1 " — + ) dt 2 + ( 1 - — + ^ ) dr 2 + r 2 dtt 2 , (2) 



and the scalar field with mass m satisfies the wave equation 

(□-m 2 )0 = O. (3) 
Resolving the field into spherical harmonics 

= £^W^)' (4) 

l,m 

we obtain a wave equation for each multiple moment; 

ttt ~ tr.r. + Vi> = 0, (5) 

where r* is the tortoise coordinate defined by 

dr * = , 02 , (6) 



1 _ 2A£ | 

r 



and the potential V is 



1(1 + 1) 2M 2Q 2 



+ 



B. The black-hole Green's function 



(7) 



The time evolution of ip is given by 

if>(r„ t)= J [G(n, <; t)Mr', 0) + G t (n, <; t)^«, 0)] rf< (8) 
for t > 0, where G(r^,r^;t) is the (retarded) Green's function satisfying 



with the initial condition G(r*,r^;t) = for t < 0. We calculate the Green's function 
through the Fourier transform; 

-r roo+ic 

G(n, <;*) = -— / G{u, <; u)e-^du, (10) 

^ J -oo+ic 

where c is some positive constant. The usual procedure is to close the contour of integration 
into the lower half of the complex frequency plane shown in Fig. [|. Then, the late-time tail 
behaviors which are our main concern should be given by the integral along the branch cut 
placed along the interval — m < uj < m. 

The Fourier component G(r*, r^; uj) in the range — m < uj < m can be expressed in terms 
of two linearly independent solutions ip\ and 1^2 for the homogeneous equation 

(£-+u?-v)$ t = 0, i = l,2. (11) 



\dr 2 , 

The boundary condition for the basic solution ipi is to describe purely ingoing waves crossing 
the event horizon, i.e., 

^~e-— , (12) 
as r* — > —00, while ip2 is required to damp exponentially at spatial infinity, i.e., 

^-e"^, (13) 



as r* — > 00, where w = \/m 2 — uj 2 . Because the complex conjugate is also a solution for 
Eq. (|TT1) , ip2 can be written by the linear superposition 



i> 2 = a^ 1 + (3^* 1 , (14) 
and the Wronskian is estimated to be 

W(uj) = V>i$ 2 ,r. - A,rM = 2iujf3. (15) 

Using these two solutions, the Green's function can be written by 



, MK,uj)ij2(n,uj) , r' > r, 

G(r*,r m ]w) = -— „ { _ _ (16) 

$i(r*,cj)V> 2 (r*,o>) , r'<r. 



2iuj(3 



The contribution G c from the branch cut to the Green's function is reduced to 



G c ( n ,rUt) = -— / Vi(r»Vi(r-*,o;)e— *dw. (17) 
Then the main task to evaluate G° is to derive the coefficients a and /?. 



III. NEARLY EXTREME LIMIT 

It is difficult in general to obtain the coefficients a and /3, since exact solutions for the 
wave equation (|TT|) cannot be expressed by any elementary function or any trancedental 
function already known. Fortunately, in nearly extreme Reissner-Nordstrom geometry, the 
procedure of asymptotic matching turn out to be very useful. Let us change the variable r 
to z defined as 



z — 1 r — r , 



2 2nr\ ' 



(18) 



where r + and r_ are the outer and inner horizon radii, respectively, and k is the surface 
gravity defined as k = (r+ — r_)/2rj_. Then we can rewrite the wave equation ( [TT]) into 



y2 ^(i 2 ^ + {2/tr 2 + z (r + + r_)} ^ 
dz 2 r +{{ z — \)kt + + 1} gLz 



-m r + {(z — l)/tr + + 1}' 



IV* {(Z - l)KT + + iy 

K 2 Z 2 -\ K+ ' 

z(r , +r H 

^ = 0. (19) 



, „ nr 2 + {2nr 2 + + z(r + + r_)} 



r 2 {(,2- l)«r+ + l} 2 
For the nearly extreme case such that 

kM < 1, (20) 
we can allows us to derive the approximate solutions valid in the region z ^> 1 or z < 1/ kM. 

A. Solutions for 2 <C 1/kM 



Expanding the wave equation as a power series in kMz and truncating terms of 
order kMz and higher, we obtain 
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+ 



dz 2 dz 
1 cj 2 1 {Au 2 M 



+ _J_ ( _ i 2u j 2 mA + {6^ 2 M 2 - 1(1 + 1) - m 2 M 2 } 
— 1 r z + 1 \ /« / 1 J 



^ = 0. (21) 



Then the solution ipi satisfying the boundary condition (|I2] ) can be written using the hy- 
pergeometric function F; 

xF ^-^ + 2^M - ^ + i -2^M - + ^ ~ + 1; ( 22 ) 
- 2 ^M- M+ i^ r(f + l)r(-2/i) 



T(f - 2iwM - /i + |)r(2io;M ~n+\] 



xF [2iuoM + //+-, — - 2iuM + \i + -, 2/i + 1; ^ 
V 2 k J 4 



to + i)r(2u) 



T(f - 2icuM + /i + |)r(2^M + fi + 1) 

xF ( 2iuM -/!+-, — - 2iuM -//+-, -2// + 1; - ) , (23) 
V 2 k 2 £/ 



where new variable £ is defined as 



(24) 



and /i is 



/i = J + + m 2 M 2 - 6cj 2 M 2 . (25) 

We have used the linear transformation formulas 15.3.6 of |J in the second equality of (p2|) . 
Using asymptotic expansions 15.7.2 and 15.7.3 of || for \iuj/k\ ^> 1, we can reduce ( j2"3"|) to 

e -^|M, (26) 

which is valid in the region 1 <C z <C 1/ kM, and will be used for asymptotic matching with 
the solutions given in z ^> 1. 
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B. Solutions for z 1 



Expanding the wave equation (|TT|) as a power series in 1/z and truncating terms of order 
1/z and higher, we obtain 

d 2 u 2 <iw 

dx 2 x dx 



+ 



W 2 M 2 - m 2 M 2 + -(4^M 2 - 2m 2 M 2 ) + \ {Qu 2 M 2 - 1(1 + 1) - m 2 M 2 } 
x x 2 



Auj 2 M 2 u 2 M 2 

+ S— + 7— + 



x 4 x(x + l) 2 
introducing the new variable x defined as 



u = 0, 



(27) 



(28) 



and the function w defined as 



ip = (x + l)u. 



(29) 



We can give the approximate solutions for Eq. (27) using trancedental functions in the 
regions i <C 1 and x 3> 1 respectively, although it is difficult to find the exact solutions 
valid in the whole range of x. 



1. Solutions for x <C 1 

Our strategy to find the approximate solutions for Eq. ( PTQ is to truncate terms of order 
and higher in the coefficients of u, since these can be smaller than other terms for i C 1, 
and to change the variable x to 



2iuM 



x 



(30) 



Then, Eq. 



can be approximated by 



d 2 u 2 du 
doc oc doc 



— A u?M 2 + \^ 2 M 2 + — \6lu 2 M 2 - 1(1 + 1) - m 2 M 2 } 

iAj %Xj JU 



u 



0, (31) 



and we can describe the solution of Eq. ([H]) using the Whittaker's functions; 

v -) 

\ x J 



/2iuM\ ( 2icuM 

u = a x W aitlJL + a 2 W- ai ^ 



x 



(32) 



where 



The asymptotic expansions for 



(7i = -2iuM. 
2iuM 



(33) 



x 



3> 1 lead to 



__. (2iuM\ _^ f.2u 



-2iuiM 



(34) 



which is necessary for asymptotic matching in the overlap region with the solutions in the 
region z <C 1/kM. On the other hand, using Eqs. 13.1.3, 13.1.4 and 13.1.33 of ||, we can 
reduce Eq. (|32|) to 

(2iuM\ 



"1 iA* 



V x ) r(|-^t-(7i 



r(-2//) ^2iuMy + ^ T(2/i) (2iuM\~^ 



T(\ + jl — <7 X ) \ x J 



(35) 



if the asymptotic expansions are applied in the region 



2 

2iuM 



x 



C 1 as an analytic extension. 



2. Solutions for x 1 



Contrary to the region i C 1, truncating terms of order x 3 and higher in the coefficients 
of u in Eq. (p7|), since these can be smaller than other terms for x ^> 1, we can approximate 
Eq. (0) by 



d 2 u 2 du 
dx 2 x dx 



+ 



u 2 M 2 - m 2 M 2 + -{Auo 2 M 2 - 2m 2 M 2 ) + -I \Qlu 2 M 2 - 1(1 + 1) - m 2 M 2 } 



u = 0. (36) 



Introducing the function Z defined as 



u 



x 



(37) 



we can reduce Eq. (^6[) into 



dx 2 



+ 



u 2 M 2 - m 2 M 2 + -(Auj 2 M 2 - 2m 2 M 2 ) + \ {Qoo 2 M 2 - 1(1 + V 



ni 



2 M 2 } 



Then we can write solutions for Eq. (|38D using Whittaker's functions; 

, M a2 ^(2wMx) M a2: ^(2wMx) 
u = bi - h b 2 — 



x 



x 



Z = 0. (3* 



(39) 



where 



a 2 = -2wM + 



m 2 M 



if estimated in the extended region \2wMx\ < 1, we obtain 

M a2 ^(2wMx) ~ (2wMxY + ^ 



(40) 



(41) 



using Eqs. 13.1.4 and 13.1.32 of ||. The solution ip 2 satisfying the boundary condition ([13] 



is 



*p2 = W a2 ^(2wMx) ~ e- mMx (2zuMx) 



\<J2 



(42) 



for \2wMx\ — > oo. 



C. Matching 



We can match both of asymptotic behaviors fl26|) at z <C 1/kM and ( |34]) at x 1 in 
the overlap region 1 C z < 1/(kM), in order to determine the coefficients ai and a 2 . On 
the other hand, we find that both of asymptotic expressions fl3"5D and (^Tj), that are the 
results due to analytic extensions from one region into the other region each other, have 
similar forms. So we can match these smoothly in order to determine the coefficients b\ and 
b 2 . In addition, considering the boundary conditions and imposed on ip% and iJj 2 , 



respectively, we can determine the coefficients a and (3 in Eq. ( |14) as follows 



a(|u;|, w) 



(3(e in \uj\,zu) 



r(2/i)r(2^+ l)(2wM)- 



r(i + /i - ( r 2 )(2|cu|M)^+5r(i + /i + 2i|w|M) 

T(-2/i)r(-2 / u + l)(2wM)^e^-5) 
r(i - /x - or 2 )(2|cj|M)-^+5r(i - fi + 2i\u\M) 



2i\u)\M 



(43) 
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and 



(3{\uj\,w) = cn(e l7r |u;|, w) 

r(2//)r(2/i+ \)(2wM)-^e^+^ 
r(| + fi - a 2 ){2\u\MY +1 2T{\ + fi- 2i\u\M) 

r(-2^)r(-2^x+ l)(2wM)^ + ie^(-^ + 5) 
Va-fi- a 2 )(2\u\M)-^ra -fi- 2i\u\M) 



-2i\u\M 

— ^ e^l*^. (44) 



IV. INTERMEDIATE LATE-TIME TAILS 

We consider the late-time behaviors of G c at the timescale 

rat > 1, (45) 

when the decaying tails will dominate. Hod and Piran |3J pointed out that for the scalar 
field with small mass, namely mM <C 1, the dominant behavior is given by Eq. ([!]) at the 
intermediate late times in the range 

mM « rat « J -L ? . (46) 

In this section we check the validity of Eqs. (^) and (0), by deriving the intermediate tail 
behaviors. Following Hod and Piran [Q, the effective contribution to the integral in Eq. ( |TTD 
is claimed to be limited to the range \u — m\ = 0(l/t) or equivalently w = 0(y/m/t). This 
is due to the rapidly oscillating term e~ lu)t which leads to a mutual cancellation between the 
positive and the negative parts of the integrand. Then, in the time scale given by Eq. (J46f) 



we note that the frequency range w = O(Wmft) leads to the inequality 

a 2 < 1. (47) 

The factor er 2 , including field's parameter ra coupled with spacetime parameter M, originates 
from the terms of order x _1 in the coefficients of u in Eqs. (|27| ) or (pB]). If the relation ( [47]) 
is satisfied, wave equation at far region can be approximated by that of flat spacetime, 
which means that the effects of backscattering due to spacetime curvature has not dominate 
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yet. In other words, the value of 02 which gives effective contributions to the integral (|PT|) 
represents a degree of the domination of the backscattering. 

The relation ( f4"T|) allows us to obtain the approximations of a and (3 as follows, 

T(2l + 1 + 2e M )r(2/ + 2 + 2t fl ) 



a(|u;|, w) 



(2wM)- l -^(2\u;\M) 



T(l + 1 + e„ - a 2 )V{l + 1 + e M + 2i\cu\M) 



-Z-l-e M _i§ (-l-l-e/j.) 



T(-2Z - 1 - 2e /1 )r(-2Z - 2e„) 
r(-Z - e M - a 2 )r(-/ - e M + 2i|w|M) 



(2roM)' +1+ ^(2|cj|M)^e^ (/+ ^ } 



x 



Hi 



2i\w\M 



-\ui\M 



(48 



and 



/3(|w|,G7) 



T(2/ + 1 + 2e M )r(2/ + 2 + 2e M ) 
r(Z + 1 + e M - a 2 )r(/ + 1 + e M - 2i|w|M) 



(2u7M)-^(2|w|M)- i - 1 - e "e^ (m+ ^ ) 



r(-2f - 1 - 2e /i )r(-2f - 2e M ) 
rH - e M - ^)r(-Z - e„ - 2i|w|M) 



{2vjM) l+1+ ^{2\uo\M) l+ ^e i ^- l -^ ) 



to 



-2i\w\M 



■MM 



(49) 



where 

e M = a* - (l + ^ O(mM) < 1. 

Expanding the ratio as a power series in mM, we can apploximate it as follows; 

a(|a;|,ct7) oc(\u\, e m w) Z! 4 



(50) 



P(\uj\,w) (3(\uj\,e™ui) (2Z)! 2 (2Z + 1)! 2 
Substituting Eq. ( |5"TD into Eq. (|X7|) , we obtain 

/! 4 



(2M) 4 '+V| 2 ' +1 ro 2Z+1 2z. 



(51) 



G c (r„r: ; t) 



(2Z)! 2 (2Z + 1)! 
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(2M) 4 ' +2 / uj 2l+1 w 2l+1 Mr*,uj)MK^)e~ iu}t du; 



+ (complex conjugate), 



(52) 



which is similar to Eq. (29) in Q, giving the damping exponent in Eq. ([!]). 
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Different from H, our analytical calculation is not based on the flat space approximation. 



The intermediate tails dominate in the range when the integrand can be approximated 
by Eq. ([47]). It is easy to find that the larger the field's mass is, the sooner it leaves the 
intermediate tails, and the phase does not appear in the case of mM > 1. 



V. ASYMPTOTIC LATE-TIME TAILS 



It is obvious from our calculation that the intermediate tail is not a final pattern of decay 
but should be replaced by another one, because the dominant contribution to the integral 
(|I~7|) is out of the region (f47|) after the intermediate late times ([46]). The change into another 
phase was also numerically suggested in [|J . Physically the change of the tail behavior will be 
a result of dominant backscattering due to spacetime curvature, which is the effect beyond 
the flat space approximation. What kind of wave pattern dominates at very late times? In 
addition, we must reveal late-time tails in the mM > 1 case of large field mass, for which 
the intermediate tails do not appear. In this section we study a tail behavior dominant at 
asymptotic late times 

1 



m z M^ 

when the effective contribution to the integral (|I~7D arises from the region 

m 2 M 



(53) 



<?2 



»1, 



W 



(54) 



different from the inequality (^) at intermediate late times. The coefficients a(tv, e ±l7T tu) 
and (3(lj, e m) are approximated for the inequality fl5"4]) by 



a(\u\,e' i7T w) = P(e iw \u\,e i7T w) 



■(2m 2 M 2 



•-2vjM , 



2vjM) 



r(2/x)r(2/x + l)(2m 2 M 2 f(2HM)-^e^-§) 
r(| + n + 2i\u\M) 

r(-2//)r(-2// + l)(2m 2 M 2 )-^(2|w|M)^-5e i i(-^-f) 



H + 2i\u\M) 



K 



2i\u\M 



\M 
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and 



(3{\u\,e- i7r w) = a{e in \u\,e i7r w) 

1 m 2 M 



e^(2m 2 M 2 ^ 



T(2fi)T(2fi + l)(2m 2 M 2 Y(2\Lu\M)-^-h^^-^ 



r(| + /i-2z|o;|M) 



r(-2/i)r(-2/i + l)(2m 2 M 2 )^(2|cj|M ) "~ s e 
r(i-/x-2i|a;|Af) 



M _i if( _ 3At i) 



3 -7r|u;|M 



Thus the contribution from this part corresponding to the dotted lines in Fig. |l| to the 
Green's function given by 



1 

Am' 



1 ' ,U ' , %(r„4(f' (|W ) e - M <i U 



(57) 



'dotted lines ^(uJjW) 

is 0(t _1 ) at most, because a(uj, e ±m w) and e ±l7r cc) converge in the limit \uj\ — > m. On 
the other hand, a(u;,tz7) and (3{uj,w) are reduced to 

(58) 

|u/|,'u;j ct^e"" |tw | , 'cu j 'fr r i/e 
for the limit ([54]), where 



a{\ui 


,zu) _ (3(e™ 




,w) _ 7 * e - i7rff2 +r/*e i7rCT2 




, w) a(e t7T 


UJ 


, w) ~~ 7e i7rcr2 + rie-i™ 2 
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r(-2/x)r(-2^ + l)(2m 2 M 2 )^(2|o;|M)^i 



r(| - fi - 2i\uj\M) 



r(2/x)r(2/x+ l)(2m 2 M 2 )^(2|^|M)-^5 , fHH .i, 



r(i + //-2z|o;|M) 



(59) 



and 



r(-2//)r(-2/i + l)(2m 2 M 2 )^(2|cj|M)^^ 



r(| - n - 2i\uj\M) 



e 2 



(-3^+4) 



r(2 / x)r(2 / i + l)(2m 2 M 2 r(2|^|M)^-l c ^ (3M+ , ) 
r(i + /i-2z|o;|M) 



(60) 



Therefore we find the contribution from this part corresponding to the dashed line in Fig. 
|l] to the Green's function to be approximated by 
1 [ m 1 a(uj,w) 



, -e l ^*V>i( r *> ^)^i( r * ) oj)duj + (complex conjugate) 

Am J uj p{uj, uj) 

Amnt 



/ e^^-^e^du + (complex conjugate), 
Jo 



(61) 
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where the phase tp is defined by 



e i V = V +7 e 

^ _|_ ^,g2l7T(T2 



and it remains in the range < </? < 27r, even if a 2 becomes very large, since we have 

fop _ | 7 |2 = _^ e ^l«|w > o. (63) 

At very late times when fl45f) and (54) are satisfied, both terms of e wt and e 2l7T(72 are 
rapidly oscillating. In physical meaning, scalar waves are mixed states with multiple phases 
backscattered by spacetime curvature, and most of these waves are canceled out each other 
by those of the inverse phase. If the value of 2na2 — tot in (|BTD is stationary at to — uq, i.e., 

— (2na 2 - tot) = 0, (64) 

particular waves with the frequency ujq remain without cancellation, and contribute domi- 



nantly to the tail behaviors. In such a case we can evaluate the integral (|61|) as the effective 
contribution from the immediate vicinity of the saddle point u . This method which is 
called the saddle-point integration allows us to evaluate accurately the asymptotic behav- 
iors of Bessel functions as a well-known example. We can find a solution for Eq. ( [64] ) 

i 



/— ~ f2irm 3 M\ 3 

w = ^Jm 2 -^ ~ I J (65) 

in the limit ujq ~ m. In order for the saddle point (|55|) to exist in the region (pl[), where 
e 2nra 2 are rapj^jy oscillating as a function of u, we need the additional relation 

mt > mM. (66) 

Approximating the integration ([H]) by the contribution from the immediate vicinity of co>o, 
we obtain 

i f i^^^y_^^^ )du ^ _^ (27r) ! (mM) i (mt) -i e - te ^o) ) (67) 



I LO = UJ : 



4irim J 4^ 
through the formula 



/oo /-oo 
cos(x 2 )dx = / sin(x 2 )dx = \ —■ (68) 
■oo J -oo V 2 
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Taking the decay rate into account, we can neglect the contribution from (|57 ) to G c in 
comparison with that from fl6"T|). Therefore tails such as ( |57|) dominate at asymptotic late 
times. Finally we arrive at the asymptotic late-time tail as 

G C {r*,r'^t) ~ -^-=(2tx)^ (mM)^ (mt)~^ sin(mt)'0i(r*, m)i/)i(r^, m), (69) 
2v 3 

of which the decay rate is independent of the location r*. 

It is found that the time, when the conditions for the application of the saddle-point 
integration, namely, (fEJ), (|53|) and (|66"D are all satisfied, must come sooner or later indepen- 
dent of field mass, and the tail (|69|) is the asymptotic behavior at late times in the limit 
rat — > oo. 



VI. SUMMARY AND DISCUSSIONS 

In this paper we have investigated mass-induced behaviors which appear in late-time 
tails of classical massive scalar fields in nearly extreme Reissner-Nordstrom background. If 
the field mass is small, namely raM <C 1, the intermediate tails given by Eq. (JI]) have been 
shown to dominate at the intermediate late-time raM <C rat <C l/(mM) 2 , consistently with 
|| (see also [||]). Our main result is the asymptotic tail with the decay rate of t~ 5 / 6 , which 
is interestingly independent of the field mass ra and the angular momentum parameter I. 
This behaviors of inverse power- law decay supports the numerical results in ||. 

Late-time tail behaviors are generally caused by the domination of the backscattering 
from far regions. It is found the asymptotic tail of massive scalar fields (^) appears when 
the effective contribution to the integral (|D]) arises from the region fl54|), namely, when the 
backscattering due to spacetime curvature dominates. Further, the frequencies of waves 
which contribute to the backscattering are sharply peaked about u . These facts suggest 
that the asymptotic t _5//6 tail is caused by a resonance backscattering due to spacetime 
curvature. 

We can also clarify the resonant picture from a viewpoint of the timescale when the t~ 5 ^ 6 
tail dominates. The basic condition for the tail dominance is given by Eq. (|45|) . However, for 
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mM <C 1, the t~ 5//6 tail requires the additional condition (|54"D , which gives mt ^> l/m 2 M 2 . 
On the other hand, for large field mass mM 3> 1, the larger mM is, the later the t -5 / 6 tail 
begins to dominate, because the tail requires the further condition ©o/m C 1 for Eq. ( |65|) 
in addition to the conditions (|45| ) and ([54]), which gives mt 3> mM. Therefore the timescale 
when the t _5//6 tail dominates will become minimum at mM ~ 0(1), which means that the 
most effective backscattering occurs for such massive scalar fields at late times. Here we 
note that the timescale when the t -5 / 6 tail dominates is determined by the black-hole radius 
M. It is conjectured that the slow decay of t~ 5//6 originates from the existence of resonant 
enhancement of massive scalar waves with the peculiar frequency u>o near the horizon. In 
conclusion, we claim that resonance behaviors in field-mass dependence are not peculiar to 
quantum scalar fields but manifest also in classical scalar fields. 

In this paper we have calculated the tail behaviors in nearly extremal limit, because we 
are motivated by the previous work . The extension of our calculation into non-extremal 
case remains in future works, to investigate whether our result in this paper is a special 
feature of the extremal case or not. 
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FIGURES 
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FIG. 1. Integration contours in the complex frequency plane. The original integration contour 
for the Green's function lies above the real ferquency axis. We choose the value of w on the dashed 
line to be w = \w\ and that on the dotted lines to be w = e ±l7T \w\. The poles in G(r^,r^;a;) are 
also shown, which give the quasinormal modes. 
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